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^. ABSTRACT 

I Based on the long-term radial- velocity surveys carried out with the McDonald Observatory 2.7ni 

C"*^ , Harlan J. Smith Telescope from 1988 to the present, we derive upper limits to long-period giant 

planet companions for 31 nearby stars. Data from three phases of the McDonald Observatory 2.7m 
planet-search program have been merged together and for 17 objects, data from the pioneering Canada- 
France-Hawaii Telescope (CFHT) radial- velocity program have also been included in the companion- 
limits determination. For those 17 objects, the baseline of observations is in excess of 23 years, enabling 
the detection or exclusion of giant planets in orbits beyond 8 AU. We also consider the possibility of 
,_!. . eccentric orbits in our computations. At an orbital separation of 5.2 AU, we can exclude on average 

O ; planets of M sin i >2.0±1.1 Mj^^p (e = 0) and M sin i >4.0±2.8 Mjup (e = 0.6) for 25 of the 31 

\^ , stars in this survey. However, we are not yet able to rule out "true Jupiters," i.e. planets of M sin 

f^ ■ i ~ 1 Mjup in 5.2 AU orbits. These limits are of interest for the Space Interferometry Mission (SIM), 

"^ ' Terrestrial Planet Finder (TPF) , and Darwin missions which will search for terrestrial planets orbiting 

"^^ nearby stars, many of which are included in this work. 

Subject headings: stars: planetary systems - extrasolar planets - techniques: radial velocities 
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C/3 . 1. INTRODUCTION What implications would the lack of such planets have 

^. ■ In the ten years since the discovery of the on theories of planet formation? How many planetary 

first extrasolar p lanet orbiting a main-sequence star systenis reseinble our own Solar system? 

llMavorfcOnelo^ lT995V more than 170 planets have In this work, we merge radial-velocity data from three 

. . . been discovered. The vast majority of these have been P^^^es of observations at McDonald Observatory, and we 

^ . detected via high-precision radial-velocity measurements ?^^^",^^ t^"" ^™™ f ^ . ^^^^ of [ Walker et alj 

■ ■ ' of the planet's gravitational influence on its host star. "IM)- In §2, we describe the data sets §3 describes 

Current Doppler surveys now routinely achieve preci- ^he merging procedure and outlines the algorithm used 

sions of 2-3 m s"! (e.g. Cochran et al. 2004), facilitat- t° ^^^^^ compamon hunts, and m §4 we present and 

ing the detection of ever-low er-mas s companions, s uch as discuss the constraints these data place on the presence 

the " hot Neptunes" ( Bonfils '^t7nfeo0 5: McArthu r"^nn "^ Pl^^^ts orbiting these stars. 

I200I iSantos et al"' 2004b; Butler ct al.l l2004|l and even 2. OBSERVATIONS 

"super-Earths" (Rivera et al. 2005). 2.1. A Brief History of Radial Velocity Planet Detection 

However, there exists a relatively long record of ra- t^ ,. i , • , r , i , 
dial velocity data at somewhat lower precision (- 15- Radial-velocity searches for planets around other stars 
20 m s-i) which ought not to be ignored. Such data began well beforej^ 1995 discovery of 51 Pegb. 
now cover nearly a quarter-century fCamnbdl fc wTIk^ | Griffin fc Griffin| 11973) showed that the fundamental 
[ia79; CaiiiDbcU ct al. 1988: Walker et al. 1995), and as l™i* *« t^^ precision of classical radial vefocity tech- 
such are extremely useful in probing nearby stars for ^iq^^s was not instruments used, but rather was a result 
long-period giant planets akin to our own Jupiter (orbital ^^ the calibration processes. They pointed out that the 
perfod 11.9 years). These data sets are valuable tools in opti^al paths followed by the stellar beam and the sepa- 
the search for extrasolar analogs to our Solar System. ^'^te comparison sourcedlmninated the spectrograph op- 
We are now beginning to obtain meaningful answers to ^'^s differently. iTulJ \m^ demonstrated how zonal er- 
some of the following questions: What is the frequency yo'^^ ™ spectrograph optics contribute directly to wave- 
of long-period giant planets in the solar neighborhood? ^"""f^ (and thus vefocity) errors. In addition, thermal 

and mechanical motion of the spectrograph m the time 

Electronic address: | robw@astro.as.utexas.edu, mike@astro.as.utexas.eda):i^ i<ag®iihaQEWigan(Q:'M.aaMng the stellar and calibration ob- 
servations also contributed to the measured velocity er- 
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rors. iGriffin &: Griffin! l)1973j) proposed that these sources 
of errors could be defeated through the cahbration of the 
spectrograph by an absorption hne spectrum imposed on 
the stellar light before it entered the spectrograph. They 
suggested that the use of the very convenient absorptions 
by the Earth's atmospheric O2 band at 6300A would 
permit stellar radial velocities to be measured to a preci- 
sion p erhaps as good as 0.01 kms~^. iCamobell fc Walken 
(|197!l) first used an externally applied absorption spec- 
trum to provide a stable velocity metric for precise radial 
velocity measurement as applied to planet detection. In- 
stead of the telluric absorption spectrum, they chose to 
use a stabilized HF gas absorption cell in front of the 
coude spectrograph of the Canada France Hawaii Tele- 
scope (CFHT). The HF 3-0 R branch hues near 8700A 
gave radial velocity precision of 15ms~^ in their twelve- 
year survey of 21 bright s olar- type s tars llCampbell et al.l 
llQSatlWalker et al.ll99iTl) . Although IWalker et all (ll99fiF 
did not claim any planet detections based on the CFHT 
data alone, two of their target objects (e Eridani and 
7 Cephei) were subsequently found to have planetary 
companions when the CFHT data were combined wit h 
McDonald Observatory data l|Hatzes et al.ll200nll2003|) . 

The CFHT survey marked the start of a large number 
of diverse efforts to conduct high-precision radial velocity 
surveys. At the University of Arizona, McMillan et al. 
l) 19851 1199 0. 1994) developed a precise radial velocity 
spectrometer based on the use of a Fabry-Perot etalon 
in transmission in front of a cross-dispersed echelle spec- 
trograph. This instrument was able to achieve velocity 
preci sion of 8ms^^ for bright sources (.McMillan ct al. 
11994(1 . and led to th e discovery of the radial velocit y 
variability of K-giants l|Smith et al.ll987HCochranll988D . 
The idea of using various types of interferometers to ob- 
tain precise stellar velocity measurements has been quite 
popular , with d esigns proposed b y Cochran ct al. (1982), 
ICoTmefi Ip85l) . lAngel X /: Wodf (1997), and Q^ (|2QQ3). 
Of these designs, onlv lvan Evken et al.. (,2003,1 were suc- 
cessful in developing the design into an instrument tha t 
actually discovered an extrasolar planet f Ce et al1l2fl0fitl . 

To our knowledge, Beckers ( 1976) first used the molec- 
ular iodine B'^HfJ',j-X^I]+ band as a velocity reference 
system, in a pplication to h is measurements of the so- 
lar rotation. iKoch fc Woeh l (1984) then extended this 
use of the B-X I2 band to the measurement of the 
Doppl er shifts o f sever al solar photospheric absorption 
lines. iLibbrechU 1(19881) then extended the use of I2 as 
a velocit y metric in searching for p-mode oscillations of 
a CMi. iMarcv & B utler (1992) then implemented the 
use of the I2 velocity metric for large-scale precise radial 
velocity surveys for extrasolar planets. The use of an 
I2 cell offers the major advantage that the I2 lines can 
be used to measure the spectrometer instrumental profile 
l(Valenti et alJl99,'il) and thus measure very precise radial 
velocitv variations ((Butler et al.lll996() . An alternative is 
to use an optical-fiber fed cross dispersed echelle spec- 
trograph with simultaneous wa velength calibration. This 
was first done with ELODIE dBaranne et al.lll996() and 
later with ELODIE ((Oueloz et al.ll2000() . Tofurther im- 
prove precision with this technique, the HARPS instru- 
ment (Rupprccht ct al. 2004) is placed in a temperature- 
stabilized vacuum chamber. 



2.2. The McDonald Observatory Planetary Search 
Program 

The McDonald Observatory Planetary Search program 
comprises a large, multifaceted investigation to detect 
and characterize planetary companions to other stars in 
our galaxy. It began in 1988 as a high-precision radial 
velocity survey of bright nearby stars using the McDon- 
ald Observatory 2.7m Harlan J. Smith Telescope and 
coude spectrograph, but has expanded substantially in 
size and scope since then. Phase I of the radial veloc- 
ity planet search program used the telluric O2 lines near 
6300A as the velocity m etric, a technique suggested by 
IGriffin fc Grlffil ((T973|) . A single order of the McDon- 
ald 2.7m telescope coude echelle spectrograph (csl2) was 
isolated onto a Texas Instruments 800 x 800 CCD at 
R = 210,000. This system gave about 20ms~^ preci- 
sion on stars down to y = 6, but suffered from sys- 
tematic velocity errors, possibly due to prevailing atmo- 
spheric winds. Diurnal and seasonal variability in the 
winds introduced spurious periodic signals in the data. 
We therefore switched to a temperature stabilized I2 cell 
as the velocity metric in 1992. This eliminated the sys- 
tematic errors, and gave a routine radial velocity pre- 
cision of ~15ms~^. This precision was limited by the 
9. 6 A bandpass of the single order of the echelle grating, 
and by the poor charge-transfer and readout properties 
of the TI 800x800 CCD. To solve these problems, and to 
achieve substantially improved precision, we began Phase 
HI of the radial velocity program in July 1998, using 
the same I2 cell with the newl y installed 2dcou de cross- 
dispersed echelle spectrograph ((Tull et al.lll994(l with its 
2048x2048 Tektronix CCD. We set up the spectrograph 
to include echelle orders from 3594A to 10762A, which 
covers the Ca II H and K lines used to measure stel- 
lar activity. Wavelength coverage is complete from the 
blue end to 569lA, and there are increa singly large inter- 
order gaps from there to the red end ((Tull et al.l ll9951. 
Using the full 1200 A bandpass of the I2 absorption band 
at the R — 60, 000 focus of the 2dcoude allows routine 
internal precision of 6-9 ms~^ to be achieved. The Mc- 
Donald program has subsequently been expanded signif- 
icantly to include a collaboration with Martin Kiirster 
and colleagues on southern hemisphere radial velocity 
surve ys using the faciliti es of ESO at La Silla and the 
VLT iJKiirster et al .1200(1 lEndl et al.l2 002: Kiirster eTall 
I2003D . and a program to use the Keck 1 HIRES spectro- 
graph to search for planetary companions to stars in the 
Hyades ((Cochran et al.ll200a iPaulson et a'Lll200l I200I 
■2004albl ) and in the use of the High Resol ution Spec- 
trogr aph on the Hobby-Eberly Telescope (Coc hran et alJ 
|2004() . Observations from the McDonald planet search 
program have been used in the dis covery of four plan- 
ets: 1 6 Cyg Bb ("Cochran et al." 1997), 7 Cep Ab 
(JHatzes et al. 2 003), e Erid ani b (Hatzes ct al. 2000), and 
HD 13189b ( Hatzes et al.ll2005|l . McDonald data were 
also instrume ntal in the detec tion of the brown dwarf 
HD 137510b l(Endl et al.ll200l . 

2.3. Observational Data Presented Here 

The majority of the data used in this study were ob- 
tained by the McDonald Observatory program Phases 
I-III using the 2.7m Harlan J. Smith telescope. A list 
of the target stars is given in Table 1. For 17 of the 31 
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stars in this study, additional data were available from 
the CF HT precision radial- velocity work of I Walker et al.l 
()1995j) . Table 2 gives a summary of all four data sets, 
including the total time span, rms of each data set, and 
the c hromospheric emission ratio log R'^k ("Novc s et alJ 
Il984|) computed from the McDonald Phase III measure- 
m.ents of the Ca II S-index. The rms values listed in 
Table 2 include the internal uncertainties of 6-9 ms^^ 
and the velocity jitter inherent to each star. 

3. DATA ANALYSIS 

3.1. Merging the Data 

Since each of the four data sets consists of veloci- 
ties measured relative to an independent and arbitrary 
zero-point, it was necessary to implement a consistent 
and robust method of combining them. Before merg- 
ing, an iterative outlier-rejection routine was applied to 
each data set separately. Data points which were more 
than 3.3 X RMS from the mean were rejected. This crite- 
rion corresponds to a 99.9% confidence level for a Gaus- 
sian distribution. The mean and RMS were then re- 
computed and the outlier rejection was repeated until 
no more points were rejected. As the targets are well- 
studied constant-velocity stars, we are confident that we 
have a well-defined distribution about the mean. For the 
visual binary 70 Oph A, several velocities were system- 
atically too low due to spectral contamination caused by 
the chance alignment of both components on the slit, and 
these points were also removed. The fainter component 
(A™ = 1.8 mag) adds a Doppler-shifted second spec- 
trum, which distorts the line shapes. Since our velocity 
computation method assumes a single set of lines, the 
resultant velocity is skewed by this contamination. 

The merging of these data sets was accomplished in 
the following manner: McDonald Phases II and III are 
joined with a trial offset, and a least-squares linear fit of 
a trend of velocity with time is performed on the com- 
bination. The offset between the Phase II and III data 
sets which minimizes the RMS about that linear fit is 
the one which is applied in order to merge two phases 
together. This process was then repeated sequentially to 
join the Phase I and then the CFHT data to the growing 
data string. The least-squares fitting allows for a linear 
trend to be present; the merging process was re-done for 
all non-binaries, this time forcing the slope to be zero 
(i.e. minimizing the RMS about the mean of the com- 
bined data for a grid of trial offsets). The RMS about 
the mean of data merged in this manner was compared 
to the RMS about a linear fit to data merged allowing 
a trend. Since our null hypothesis is that the stars are 
radial-velocity constants, the method which allowed a 
slope was only chosen for the 9 stars which showed an 
improvement in the RMS of the merged data set when 
a slope was allowed. Those stars were: 77 Cep, 77 Cas, a 
For, e UMa, ^ Boo A, /x Her, 70 Oph A, and 61 Cyg A & 
B. All of these stars are well-known long-period binaries, 
so the use of a slope was justified. 

The orbital periods of these 9 known binaries were suf- 
ficiently long that a simple linear approximation was suf- 
ficient, except for /i Her, ^ Boo A (quadratic trend) and 
70 Oph A (Keplerian solution) . Fits to these trends were 
subtracted from the data sets and the residuals were used 
as input for the detection-limit algorithm. The three 
data sets on the visual binary /i Her were merged in a 



similar manner. The binary period (43.2 years; Worley 
& Heintz 1995) is such that a linear trend is a subopti- 
mal approximation for the purposes of merging the 16.9 
years of data, and a quadratic fit was used to approx- 
imate the shape of the orbit. This procedure was also 
used for ^ Boo A; the data and fits are shown in Fig- 
ure 1. For 70 Oph A, the quadratic approximation was 
still insufficient, as the 16.1 years of available data rep- 
resents a signifi cant portion of the 91.44 year published 
orbital period ("Ba tten et alJll98^ . and an orbital so- 
lution for the binary system had to be obtained using 
GaussFit ij.leffervs et al.lll987l) . The solution to the com- 
bined data set (McDonald Phases I-III plus much lower- 
precision data from Batten et al. 1984) is given in Table 
3, and plotted in Figure 2. The large uncertainties in the 
fitted parameters are due to the fact that the available 
data do not encompass a complete orbit, and hence the 
model is poorly constrained. The fitted value of w was 
derived by a rudimentary grid search around the pub- 
lished value, as GaussFit would not converge with oj as 
a free parameter. The orbital parameters of the binary 
system are substan tially different than those reported in 
iBatten et al.ll)1984D . likely owing to the fact that the Mc- 
Donald velocities, which are an order of magnitude more 
precise, drive the fit. Only the McDonald data were used 
in the computation of companion limits for 70 Oph A. 

3.2. Periodogram Analysis 

After the data sets were merged together, and trends 
due to binary orbits were subtracted, we searched 
for periodic i ties using a Lomb-Scargle periodogram 
IjLombl IT976t [Scarglc 1982). False alarm probabihties 
(FAP) were established usin g a bootstrap rando miza- 
tion method as d escribed in iKiirster et al.l l)1997f) and 
lEndl et all l)2002() . This bootstrap method, unhke the 
nominal analytic FAP formula of Scarglc (1982), makes 
no assumptions about the distribution of the data; 
rather, the error distribution of the actual data them- 
selves is used. The periodogram search interval was from 
2 days up to the full extent of observations for each ob- 
ject. We used 10,000 bootstrap randomizations to deter- 
mine the FAP of the highest peak in the periodogram for 
each star; the results are given in Table 4. Only one star 
had a periodicity with a significance greater than 99%: 
TT^ Ori, with a peak at 73.26 days and a bootstrap FAP 
of 0.8%. As this is an active early-type (F6V) star, we 
are confident that this is due to a combination of stellar 
activity and reduced velocity precision (from the paucity 
of lines) rather than the presence of a companion. Ad- 
ditionally, no peak was evident in a periodogram on the 
more precise McDonald Phase III data alone, supporting 
our conclusion that it is a spurious noise spike. 

3.3. Determination of Companion Limits 

Companion limits were determined via an algorithm 
which injected test signals into the data, and then at- 
tempted to recover that signal using a periodogram 
search. Test signal s were genera ted using the method 
of Lehmann-Filhes llLehmann-Filheslll894l) as described 
in lBinnendiikI 1)1960(1 : 

Vr = ^ + K[e cos uj + cos{v -I- w)] , (1) 

where Vr is the radial velocity, 7 is the systemic veloc- 
ity, K is the velocity semiamplitude, e the orbital ec- 
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centricity, u the argument of periastron, and v the true 
anomaly. In the above equation, ly can be expressed in 
terms of observables via the foUowing relations: 



tan — = 
2 



(1 + e) 



(1-e) 



1/2 E 

tan — 
2 



and 



M : 



27r, 



To) = E ~ e sin E, 



(2) 



(3) 



for a signal with period P, periastron passage time To, 
observation time t, mean anomaly Af, and eccentric 
anomaly E. For each data set, the algorithm stepped 
through 300 trial periods at even steps in the logarithm 
between 2 days and the total duration of observations. 
At each trial period, the program sampled 30 phase steps 
of the periastron time Tq, and for nonzero eccentricities, 
sampled 36 values of the argument of periastron w at in- 
tervals of 10 degrees. The velocity semi-amplitude K of 
the test signal was allowed to vary from 10 to 1000 m 
a~^. The systemic velocity 7 of each combined data set 
was forced to be equal to zero. For each set of Keple- 
rian orbital parameters, synthetic radial velocities were 
generated using the observation times from the input 
data. This simulated signal was then added to the data, 
and a Lomb-Scargle periodogram was used to attempt 
to recover that signal. For a signal to count as hav- 
ing been recovered, the pcriodogram's highest peak had 
to occur at the correct period with a FAP of less than 
0.1%. The FAP was es timated using the formula from 
IHorne fc BaliunasI (1986). For each K velocity, 99% of 
the 1080 test signals^ had to be recovered in this manner 
in order for that velocity semiamplitude (corresponding 
to a planet of a given mass at a given semimajor axis) to 
be considered "ruled out" by the data. If more than 1% 
of the orbital configurations tested at a given K value 
were not ruled out in this manner, the algorithm in- 
creased K hy 1 m s^^ and repeated the process until 
it was able to recover 99% of the parameter configura- 
tions (e,a;,To) that were tested. Once this occurred, the 
planetary mass ruled out was computed by the following: 



M2 sini = (1-e. ) 



2N1/2 



(1.036 X 10-'^)MlPK" 



1/3 



(4) 



where Mi is the stellar mass in solar masses, P is the or- 
bital period in days, K is the radial-velocity semiampli- 
tude in km s~^, and AI2 sin i is the projected planetary 
mass in solar masses. The algorithm then moved to the 
next trial period and repeated the entire process. 

Unlike most previous comp a nion limit d e termi- 
nations j Murdoch et al.' "1993; "Walker et al." 119951 
Nelson fc A ngel 1998; Gumming et al. 1999; End l et all 
2002J) . this procedure allows for nonzero eccentricities 
in the trial orbits (but see Desidera et al. 2003). For 
the eccentric case, a value of e = 0.6 was c hosen; of the 
known extrasolar planets, 90% have e < O.G lMarcv et al.l 
((2P05). Allowing higher eccentricities substantially 
reduced the ability to rule out the test signals, as the 
sporadic sampling of the data would likely miss the 

^ For signals with e = 0, all 30 variations in Tq had to be recov- 
ered, as w had no meaning. 



relatively high-velocity points near periastron. Figure 3 
shows the effect of allowing various ranges of eccentricity. 
Note that although the limits derived by this study are 
rendered somewhat less stringent by the inclusion of 
nonzero eccentricities, the effect is relatively minor for 
our adopted value of e = 0.6. Allowing a larger range 
of eccentricities (up to e = 0.9), however, substantially 
reduces the sensitivity of the companion-limit determi- 
nation, as demonstrated in Figure 3. We emphasize that 
the upper limits derived by the method described above 
are much more stringent, and hence will result in higher 
companion-mass limits than those reported by previous 
studies. 

4. RESULTS AND DISCUSSION 

Limits to planetary companions derived using these 
data are shown in Figures 4-11. In each panel, the lower 
set of points (solid line) represents the companion limits 
for the zcro-cccentricity case, and the dotted line is for 
the case of e = 0.6. Notably, despite the abundance 
and quality of data available in this study, we are as yet 
unable to rule out any planets with M sin i ^ 1 Mj^p 
in 5.2 AU orbits with eccentricities as large as e = 0.6. 
When only circular orbits are considered, such objects 
can be ruled out for r Ceti, a Dra, 61 Cyg A, and 61 Cyg 
B. Table 5 lists the minimum planet masses that can be 
excluded by these data at selected semimajor axes, for 
the e — and e — 0.6 cases. The results given in Table 
5 are shown in histogram format in Figure 12, which 
indicates that for most stars in this survey, Saturn-mass 
planets in close orbits (a ~ 0.1 AU) can be ruled out. 

It is also useful to consider the effect of giant planets 
in intermediate orbits (a ~ 2 — 3 AU) which may perturb 
lower-mass planets within the habitable zone of the star. 
If such objects can be excluded with confidence, their 
host stars become attractive candidates for the Terres- 
trial Planet Finder (TPF) and Darwin missions, which 
aim to det ect Earth-like planets in the habitable zone. 
IMenou fcT abachnik (2003) defined a planet's zone of 
influence to extend from Rm = (1 — e)a — ZRhui to 
Rout = (1 + e)a + SRhui, where the Hill radius is 



Rhui 



3AL 



1/3 



(5) 



and e is the planet's eccentricity, a is its semimajor axis, 
Mp is it s mass, and M^ is the mass of the star. Simula- 
tions bv IMenou fc TabachnikI l)2003l ^ demonstrated that 
terrestrial planets were nearly always ejected or con- 
sumed in systems where an eccentric giant planet's zone 
of influence overlapped the habitable zone. We can then 
ask whether the companion limits derived in this work 
can be used to exclude such perturbing bodies. Such 
a pursuit is limited by the fact that even distant giant 
planets can disrupt the habitable zone if their orbits are 
sufficiently eccentric, and as shown in Figure 1, the na- 
ture of the radial-velocity data is such that we are least 
sensitive to the most eccentric planets. Nevertheless, for 
small eccentricities (e^O.2), it is possible to combine 
these dynamical calculations with our companion-limit 
determinations to define a "safe zone" : a region of pa- 
rameter space in which we can exclude perturbing giant 
planets exterior to the habitable zone. For such regions. 
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the possibility of terrestrial planets in the habitable zone 
remains open for programs such as TPF and Darwin. In 
Figures 4-11, the region to the left of the dot-dashed line 
and above the dotted line defines the "safe zone" for per- 
turbing outer giant planets with e < 0.2. These were only 
plotted for main-sequence stars, using th e definition of 
the "c ontinuously habitable zone" given in lKasting et al.l 
l)1993j) . The region left of the dot-dashed line and be- 
low our limits represents a set of potentially dangerous 
objects which would disrupt the habitable zone, yet be 
undetectable with the current data. Higher eccentricities 
would push the dot-dashed line to the right and reduce 
its slope, such that for perturbers with e > 0.5, our limits 
computations can say nothing about such objects (i.e. the 
curves would not intersect). The companion limits we 
have derived thus place some constraints on potentially 
disruptive objects in these systems, which will assist in 
target selection for the TPF and Darwin missions. 

Noting that the Phase III data are of substantially 
higher quality than the previous data sets, we asked what 
velocity rms would be required to rule out a Jupiter ana- 
log orbiting a solar-type star. We generated simulated 
observations consisting of Gaussian noise at the actual 
observation times (spanning 16 years) for 16 Cyg A, 
the star in this study which is closest in spectral type 
(G1.5 V) to our Sun. Figure 13 shows the results of 
the companion-limit algorithm on four of these simulated 
data sets with four levels of rms scatter. In order to rule 
out a planet with M sin z of 1 Mjup in a 5.2 AU or- 
bit (e = 0.1), the data need to have an rms less than 
about 10 m s~^. The rms of Phase III observations of 
16 Cyg A is 6.9 m s~^ over a period of 6.3 years, whereas 
the complete 16 years has an overall rms of 26.8 m s~^. 
Hence, a Jupiter analog could be ruled out for 16 Cyg A 
with about 10 more years of data of the same quality as 
McDonald Phase HI. The weighted mean rms of all Mc- 
Donald Phase HI observations in this survey is 12.6 m 
s~^. Eleven stars (see Table 2) currently achieve a Phase 
HI rms better than 8 m s~^; this represents one-third 



of the stars discussed in this work. These simulations 
show that the average precision of Phase HI needs to be 
improved by about 2-3 m s~^ in order to achieve the sen- 
sitivity required to detect or exclude Jupiter analogs for 
all of these stars. 



5. SUMMARY 

We have amassed a substantial quantity of radial- 
velocity data on 31 nearby stars, spanning up to 23 years 
of observations. We have applied a robust method to 
place conservative limits on planetary companions orbit- 
ing these stars. We have considered the effect of eccen- 
tric orbits in our computations, which better reflects the 
diversity of known extrasolar planetary systems. Saturn- 
mass planets within 0.1 AU can be ruled out for nearly 
all of these stars. At 5.2 AU, we can exclude on average 
planets with M sin i >2.0±1.1 Mjup (e = 0)and M sin 

i >4.0±2.8 Mjup (e = 0.6) for 25 of the 31 stars. Al- 
though we now have a quite sufficient time baseline for 
the detection of Jupiter analogs in ~12 year orbits, we 
find that the overall velocity precision is not yet sufficient 
to exclude Jupiter-mass planets at 5.2 AU. However, im- 
provements in the current McDonald Observatory 2.7m 
telescope planet search program put the desired level of 
precision within reach for inactive stars. 
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TABLE 1 
Target List and Stellar Parameters 



Star 



HR 



Spec. Type 



V magnitude 



Mass (M0) 



Reference for Mass Estimate 



r] Cas 
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GOV 


3.44 


T Get 
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G8V 


3.50 


e Per 
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F8V 
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i Per 
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GOV 


4.04 


a For 


963 


F8V 
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«! Get 


996 


G5V 


4.82 


5Eri 


1136 


KOIV 


3.54 


o^Eri 


1325 


KIV 


4.43 


7r3 Ori 


1543 
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3.19 


A Aur 
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e UMa 
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3.17 
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/3 Vir 


4540 


F8V 


3.61 


/3 Com 
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A Ser 
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7 Ser 
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36 Oph A 
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5.29 


fj, Her 


6623 


G5IV 


3.42 
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a Dra 
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72 Del 
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L PSC 
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4.13 
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TABLE 2 
Summary of Observations 



Star 


N 


T 


GFHT rms 


Phase I rms 


Phase II rms 


Phase III rms 


1°S R'hk 






(years) 


(rns-i) 


(m s-i) 


(m s-i) 


(m s-l) 




■q Gas 


131 


16.3 




28.1 


28.9 


7.3 


-4.926 


rCet 


183 


23.3 


14.2 


22.3 


28.7 


10.4 


-4.979 


Per 


65 


9.5 






59.6 


17.9 


-4.919 


L Per 


165 


23.9 


18.2 


30.0 


18.6 


10.2 


-5.041 


a For 


66 


15.3 




28.9 


41.3 


15.8 


-5.023 


«! Get 


134 


23.0 


23.7 


34.4 


29.7 


20.1 


-4.441 


5Eri 


109 


16.1 




17.8 


22.5 


9.4 


-5.228 


o2 Eri 


138 


22.1 


18.6 


30.4 


18.0 


14.6 


-4.951 


7r3 Ori 


160 


15.9 




169.9 


113.3 


26.4 


-4.716 


A Aur 


63 


10.4 






20.6 


10.1 


-5.051 


e UMa 


268 


18.3 


18.8 


43.1 


43.7 


14.8 


-5.608 


36 UMa 


194 


23.0 


21.0 


19.4 


22.4 


8.2 


-4.811 


/3 Vir 


205 


23.3 


28.4 


30.5 


32.6 


7.6 


-4.942 


/3 Gom 


191 


22.4 


18.4 


31.7 


42.0 


11.1 


-4.749 


61 Vir 


149 


23.0 


18.4 


25.8 


31.0 


7.8 


-5.030 


^ Boo A 


186 


23.3 


23.6 


34.5 


31.9 


23.3 


-4.420 


A Ser 


63 


11.5 






7.9 


19.4 


-4.936 


7 Ser 


149 


15.1 




84.7 


41.6 


25.6 


-4.934 


36 Oph A 


91 


21.9 


20.1 


21.5 


33.1 


16.3 


-4.614 


/I Her 


173 


16.9 




28.4 


24.4 


9.2 


-5.092 


70 Oph A 


98 


16.2 




111.4 


43.3 


17.4 


-4.736 


a Dra 


178 


23.6 


14.5 


21.5 


23.1 


9.9 


-4.865 


16 Gyg A 


102 


16.0 




34.4 


29.0 


6.9 


-5.018 


31 Aql 


38 


7.9 






22.5 


11.9 


-5.123 


/3AqI 


183 


22.1 


14.6 


28.0 


20.1 


16.4 


-5.171 


7^ Del 


103 


15.0 




23.7 


21.7 


16.8 


-5.354 


r; Gep 


187 


23.5 


19.2 


29.5 


16.4 


9.6 


-5.223 


61 Gyg A 


143 


23.4 


20.7 


22.3 


13.5 


7.2 


-4.862 


61 Gyg B 


121 


22.4 


16.9 


23.5 


16.2 


3.9 


-4.962 


HR 8832 


119 


22.9 


14.9 


22.8 


13.8 


12.9 


-5.013 


t Psc 


54 


11.2 






26.4 


10.3 


-4.915 



TABLE 3 
Orbital Solution for 70 Oph A 



Parameter 



Estimate 



Uncertainty 



Period (years) 

To (HJD) 

e 

LU (degrees) 

Ka (m s-l) 



179.1 

2445610.22 

0.69 

193.6 

4163 



99.5 

91.74 

0.11 

257 
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TABLE 4 
Results of Periodogram Analysis 



Star Period (days) FAP 

T] Cas 2.964 0.296 

T Get 8.389 0.034 

Per 348.432 0.443 

L Per 91.996 0.027 

a For 8.262 0.395 

K^ Get 33.069 0.034 

5 Eri 22.051 0.788 

o2 Eri 7.463 0.900 

■K^ Ori 73.260 0.008 

A Aur 2.549 0.188 

9 UMa 2.587 0.178 

36 UMa 3703.704 0.694 

13 Vir 3.307 0.262 

13 Com 9.946 0.315 

61 Vir 9.777 0.072 

S, Boo A 6.298 0.202 

A Ser 4.929 0.773 

7 Ser 8.961 0.137 

36 Oph A 2.641 0.565 

11 Her 7.623 0.039 

70 Oph A 6.138 0.989 

a Dra 30.553 0.109 

16 Gyg A 2.873 0.038 

31 Aql 8.374 0.536 

13 Aql 2.800 0.019 

7^ Del 3.837 0.099 

■q Cep 5.940 0.030 

61 Gyg A 6.570 0.159 

61 Gyg B 6.872 0.846 

HR 8832 25.227 0.144 

L Psc 2.499 0.144 
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TABLE 5 

Minimum-Mass Companion Limits 



Star 


M sin i (Mjup) 


M sin i (Mjup) 


M sin i (Mjup) 


M sin i (Mjup) 


M sin i (Mjup) 


M sin i (Mjup) 




0.05 AU 


0.1 AU 


3 AU 


3 AU 


5.2 AU 


5.2 AU 




e = 0.0 


e = 0.0 


6 = 0.6 


e = 0.0 


6 = 0.6 


6 = 0.0 


ri Cas 


0.13 


0.17 


1.30 


0.87 


2.31 


1.37 


rCet 


0.09 


0.16 


1.14 


0.67 


1.34 


0.90 


e Per 


0.48 


0.62 


5.71 


3.46 






t Per 


0.12 


0.16 


1.46 


0.95 


2.70 


1.69 


a For 


0.32 


0.57 


6.12 


3.05 


8.46 


4.37 


k1 Get 


0.18 


0.27 


2.34 


1.44 


3.81 


2.05 


5Eri 


0.13 


0.19 


2.10 


1.07 


12.74 


1.88 


0^ Eri 


0.12 


0.16 


1.33 


0.88 


1.66 


1.16 


7r3 Ori 


0.84 


1.51 


9.79 


6.01 


46.70 


8.54 


A Aur 


0.18 


0.24 


2.19 


1.25 






e UMa 


0.24 


0.36 


3.09 


1.96 


4.60 


2.57 


36 UMa 


0.13 


0.17 


1.46 


0.91 


2.54 


1.71 


/3 Vir 


0.16 


0.23 


2.49 


1.57 


4.19 


2.61 


13 Com 


0.16 


0.26 


1.95 


1.24 


3.29 


2.22 


61 Vir 


0.14 


0.20 


1.62 


1.16 


2.58 


1.69 


S, Boo A 


0.16 


0.19 


2.00 


1.24 


2.61 


1.78 


A Ser 


0.16 


0.21 


2.01 


1.14 






7Ser 


0.44 


0.49 


4.15 


2.97 


10.89 


5.80 


36 Oph A 


0.19 


0.23 


2.24 


1.39 


5.83 


2.33 


/I Her 


0.14 


0.21 


1.53 


0.96 


2.55 


1.52 


70 Oph A 


0.46 


0.84 


8.71 


5.04 


12.78 


7.19 


cr Dra 


0.11 


0.12 


1.12 


0.79 


1.60 


1.03 


16 Cyg A 


0.18 


0.32 


2.21 


1.38 


5.44 


2.45 


31 Aql 


0.22 


0.38 


... a, 


1.90 






/3Aql 


0.12 


0.18 


1.61 


1.04 


2.67 


1.77 


72 Del 


0.22 


0.32 


2.55 


1.60 


4.70 


2.32 


T] Cep 


0.13 


0.17 


2.02 


0.86 


2.41 


1.52 


61 Cyg A 


0.09 


0.14 


1.60 


0.85 


2.10 


0.98 


61 Cyg B 


0.07 


0.10 


1.12 


0.66 


1.48 


0.80 


HR 8832 


0.11 


0.14 


1.65 


0.81 


2.56 


1.14 


I. Psc 


0.26 


0.38 


4.35 


1.85 


4.91 


2.64 



^Too few data points for a reliable periodogram search, due to undersampling of eccentric test signals. 
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Fig. 1. — Quadratic fits to /i Her data (left panel) and ^ Boo A (right panel). 
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Fig. 2. — Orbital solution for the visual binary 70 Oph A using data from lBatten et aLI 119841) (filled circles) and McDonald Observatory 
phase I-III (open circles). 
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Log Period (days) 



Fig. 3. — The effect of eccentricity on limit determinations. Allowing higher eccentricities reduces the sensitivity somewhat, due to the 
increased probability of unfortunately-phased observations. 
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Fig. 4. — Planetary companion limits for r) Cas, r Get, 6 Per, and t Per. The vertical dashed line indicates the 11.87 yr orbital period 
of Jupiter. Planets in the parameter space above the plotted points are excluded at the 99% confidence level. The solid line represents 
the companion limits for the zero-eccentricity case, and the dotted line is for the case of e = 0.6. The region left of the dot-dashed line 
represents the set of perturbers at e = 0.2 which would disrupt the habitable zone. 
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Fig. 5. — Same as Figure 4, for a For, k^ Get, <5 Eri, and o^ Eri. 
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Fig. 6. — Same as Figure 4, for n^ Ori, A Aur, UMa, and 36 UMa. 
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Fig. 7. — Same as Figure 4, for f3 Vir, /3 Com, 61 Vir, and g Boo A. 
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Fig. 8. — Same as Figure 4, for A Ser, 7 Ser, 36 Oph A, and /i Her. 
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Fig. 9. — Same as Figure 4, for 70 Oph A, a Dra, 16 Cyg A, and (3 Aql. 
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Fig. 10. — Same as Figure 4, for 31 Aql, 7^ Del, ij Cep, and 61 Cyg A. 
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Fig. 11.— Same as Figure 4, for 61 Cyg B, HR 8832, and t Psc. 
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Fig. 12. — Histogram showing the lowest M sin i of planets ruled out at seinirnajor axes of 0.05 AU (left panel) and 0.1 AU (right panel). 
The dotted vertical line indicates the mass of Saturn, and the dashed line indicates the mass of Neptune. 
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Fig. 13. — Planetary companion limits for 16 Cyg A, using simulated data matching our observation times, but with an rms of 14 m 
s~^ (red, open circles), 11 m s~^ (green, solid line), 8 m s~^ (blue, filled circles) and 5 m s~^ (magenta, dotted line). The vertical dashed 
line indicates the 11.87 yr orbital period of Jupiter. Planets in the parameter space above the plotted points are excluded at the 99.9% 
confidence level. This shows that a Jupiter analog could be ruled out at the 99.9% level with a sufficiently long baseline of data similar in 
quality to the McDonald Observatory phase 3. 



